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In this work, we focus on some properties of the parity-even sector of the CPT-even electrody- 
namics of the standard model extension. We analyze how the six non-birefringent terms belonging 
to this sector modify the static and stationary classical solutions of the usual Maxwell theory. We 
, observe that the parity-even terms do not couple the electric and magnetic sectors (at least in the 

stationary regime). The Green's method is used to obtain solutions for the field strengths E and 
B at first order in the Lorentz- covariance-violating parameters. Explicit solutions are attained 
for point-like and spatially extended sources, for which a dipolar expansion is achieved. Finally, 
it is presented an Earth-based experiment that can lead (in principle) to an upper bound on the 
. anisotropic coefficients as stringent as (ke-)'-' < 2.9 x 10"^". 

PACS numbers: ll.30.Cp, 12.60.-i, 41.20.-q, 41.20.Cv 
^! I. INTRODUCTION 

^, 

In recent years, investigations concerning Lorentz symmetry violation have been undertaken mainly in the 
, context of the standard model extension (SME) developed by CoUaday and Kostelecky 0, which incorpo- 
^ ' rates Lorentz-invariance-violating (LIV) terms in all sectors of the usual standard model of the fundamental 
. interactions. The abelian or electromagnetic sector of the SME is composed of a CPT-even and a CPT-odd 
QQ \ part. The CPT-odd sector is represented by the Carroll-Field- Jackiw term, epap^pV^ A"' F'"-'' , whose properties 
were first examined in Ref. [^. The investigations on this electrodynamics have been performed in a broad 
, perspective, addressing aspects as diverse as the consistency and quantization of the model radiative cor- 
^—1 ■ rections Q , classical solutions @] , Cerenkov radiation , cosmic background radiation 0] , and other features 
00 . 0. More recently, the CPT-even sector, represented by the term WaiypLpF°"^ PP'-f , has been investigated as well, 
' embracing the study of small deviations of the Maxwell electrodynamics stemmin g fr om this term and some 
attempts of imposing upper bounds on the LIV parameters 111, 12 1, [l3|, 14 1 . |l5| . 16|, 17|. 



> 
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The Lagrangian density of the CPT-even electrodynamics of the Standard Model Extension has the form 

1„ 1 



?^.\ L = --F^.F^'' - -Wc^.p^F^-^FP"^ - Jo^A^, (1) 



where the background tensor Wavp^p has the same symmetries as the Riemann tensor 
\Waypv = -W^ap^, Wayp^ = -Wau^p, Wayp^ = Wp^au] and a double null trace, W^afS = 0, implying 
19 components. This tensor Wavptp can be written in terms of four 3x3 matrices kde,kdb, khe, khb, defined 



in Refs. PJ, DJ as 



The matrices kde,i^hb contain the parity-even components and possess together eleven independent compo- 
nents, while k,db,hhe possess together eight components and describe the parity-odd sector of Wavp'^- Four 
tilde matrices and one trace element can be written as suitable combinations of kdb, kdb, i^he^ k-hb'- 

{Tie+y - li^DE + ^iHBy\ (^e-)^' - ^(^^DB ^ - ^P^^^DeT, (3) 

[Ko+Y'' = ^{kdB + HHeY^, {Ko-Y'' = \il^DB - KheY'', Ktr = ^{kdeT' ■ (4) 
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From the eleven independent components of the matrices Ke+,Ke , the five elements enclosed in Ke+ are 
constrained by birefi'ingence to the level of 1 part in 10"^^ (see Refs. 13|,[l4|), there remaining six non-birefringent 
ones (the trace element and the five components of the matrix He-) to be constrained by other methods. 
From the eight elements of the parity-odd sector, five (contained in the matrix Ho-) are tightly bounded by 
birefringence, there remaining only three components (belonging to Ko+), which were parameterized as the k 



vector 15|, written as = ^e^P"^ (i^db)^'' ■ In some recent papers 16|, the absence of Cerenkov radiation from 
ultrahigh-energy cosmic rays (UHECRs) has been used to state bounds at the level of 1 part in 10^^ on the nine 
nonbirefringent terms of W^^p^, belonging both to the parity-even and parity-odd sectors. 

In Ref. there was performed an analysis focused on the three non-birefringent components (k^) of 

the parity-odd sector of W (the parity-even components were taken as null in order to isolate the parity-odd 
sector physics) . The stationary classical solutions for the Maxwell electrodynamics modified by these three LIV 
coefficients were properly evaluated by means of the Green's method. With these solutions, it was described a 
device able to yield a nice upper bound, < 10^^^, in the context of an Earth-based experiment. 

The aim of the present work is to study the stationary aspects of the classical electrodynamics stemming from 
the parity-even sector of the tensor Wavpip- For that, we use the following parameterization kde = —khb, and 
we consider as null the parity-odd sector. The goal is to determine how the the six non-birefringent parity-even 
components modify the classical and stationary solutions for Maxwell electromagnetism. Certainly, the idea is 
also to use the results obtained to properly constrain the magnitude of the LIV coefficients. 

This work is organized as follows. In Sec. II, we write the wave equations and modified Maxwell equations 
and apply the Green method in order to obtain the required stationary solutions. In Sec. Ill, we present our 
final remarks and describe a measurement device able to yielding an upper bound on the LIV parameter as 
stringent as (^e-)'^ < 2.9 x 10"^°. 



II. WAVE EQUATIONS AND STATIONARY CLASSICAL SOLUTIONS 



Here, the focus is on the classical properties of the parity-even part of the tensor W, particularly on those 
produced by the six non-birefringent components (located in the matrix Ke_ plus the trace element Ktr)- Hence, 
we take the parity-odd sector as null {kdb = hhe = 0) to isolate the physics of the even sector. Further, we 
adopt the following parameterization ( kde ) = — (khb) , which implies (Ke+) = 0, by considering the stringent 
bound imposed by birefringence data 3, 14 1. Moreover, it implies the following relation for the non-birefringent 
components: 



{KDEy'' - n5^\ n = -6^''Tr{KDE)- (5) 
In order to evaluate the classical solutions of this model, we write the wave equation for the four-potential 

UA°- - 2W°'''P^d^dpAx = J", (6) 
which yields two differential equations, one for the scalar potential, and one for the vector potential, 

"(1 + n)d1 - (1 + n)v2 - (k,_)'^' -t- (Ke-)'"' d^dtA.j = p, (7) 

[(1 + n)dl - (1 - n)V2] A, - 2ndtd,AQ - {H^-)"' dtE, - (^e^)^'' d^Bi = j,. (8) 



where we have used Ej 
are read 



2 ^ VJ ^ P3 ' 



(Ke-),^ ■ At the stationary regime, such equations 



(1 + n)V^ + (He-)'' d,d, 



Ao 



[(1 - n)W^] A, + (Ke-)^'' dj,Bi = 



(9) 
(10) 



These equations reveal that the electric and magnetic sectors are decoupled (in the stationary regime) in 
contrast with the electrodynamics of the parity-odd sector, in which these sectors are entirely entwined (see 



3 



Ref. Applying the differential operator eaudb to Eq. (fTU]) . we obtain the following differential equation 

for the magnetic field 



((1 - n)5al - {T^e-T') + {Ke-Y' dadj] Bl ^ ~ [V X j) ^ 



(11) 



While the homogeneous Maxwell equations remain unmodified ( V x E + 9tB — 0, V-B — 0), the inhomogeneous 
ones (Gauss and Ampere law) are altered, taking the form 



(l + n)V-E- (Ke-)''a,^j =p, 



(12) 
(13) 



(1 + n)dtE, - (1 - n) (V X B),^ + e„, (^e-) d^Bi + {Ke-Y" dtEq= -j,. 
In the stationary regime, the latter equation provides 

(1 - n) (V X B)^ - e,,, (Ke-)"' d,Bi^]„ (14) 
which under the action of the operator curl operator {eaudb) yields the same expression as Eq. (jlip . 



A. The Green's function for the scalar potential 



The solution for the scalar potential may be obtained by the Green's method. The Green's function for Eq. ^ 
fulfills 



(1 + n)V^ + (Ke-),^. d,dj\ G(r - r') = ^^(r - r'), 
and the scalar potential is given as 



(15) 



(16) 



The Green's function in Fourier space is given as G{r — ?-') = {2tt) j d^p G (p) exp [— ip • (r — r')], so that 
we obtain 



1 - n - (Ke^y 



PtPj 



(17) 



at first order in the LIV parameters. Remembering that the LIV coefficients are small, we 



1 + n+ {Ke^y^ PiPj/p'^ 



used 

the Green's function takes the following form: 



1 — n — {Ke-Y'' PiPj /p^ ■ Carrying out the inverse Fourier transform, 



G(r - r') 



1 

47r 



(l-n) 



r — r 



2|r-r'| 



(18) 



It presents a genuine Coulomb contribution screened by the factor (1 — n) and a non-Coulomb contribution 
related to the LIV non-isotropic coefficients (Ke_)j^-. The overall behavior, r~^, remains the same as happens 
in the Maxwell electrodynamics. 

Using the Green function and Eq. the scalar potential due to a general charge distribution [p (r')] 

is 



^0 (r) 



1 

47r 



(l-n) / dh' + (fi^.r I d't 



, {v -r')i{r -v')j 



which implies the following electric field strength: 



(r) 



1 

47r 



(l-n) / dVp(r')^^— ^-(«e-)'^ 



2 r — r' 



(19) 



j dh'p (r') 



3 (Tie- f J d'rV(r 



, , „ (r-rO,(r-rO^.(r-rO' 



2|r-r'|^ 



(20) 



4 



With this expression, we may immediately evaluate the scalar potential and electric field strength for a 
point-like charge at rest [p{r') = qS{r')], yielding 



^0 (r) 
E'^ (r) 



_<1_ 

An 



(21) 
(22) 



The scalar potential and the electric field present a genuine Coulomb contribution, with the screening factor 
(1 — n) , and a non-Coulomb contribution related to the LIV coefficient (ke-)*"' • This latter term leads to 
variations of the scalar potential and electric field along a circular path around the point-like charge. Such an 
effect can be used to impose an upper bound on the LIV parameters, as will be described in Sec. III. 

From the expressions (I19|) and (j20p . for an arbitrary charge distribution, we express the scalar potential and 
the electric field in the dipolar approximation (|r — r'|~^ = + r • r' /r^): 



E^«-:^^(i-«) 



H 5-^ 



2r3 



PI 3(r.Pe) 



+ 3(Ke-) ,5 



I3 ^ 

2r5 



qr^ - PI + 



r Pe 

rirj(q + 3 — j^) - nPej - rjPe, 



(~ VI ^ ( P ^ o (l--Pe) 



3 [Ke-) 



(23) 



(24) 



where q = J p{v')d?'v' and Pe=: / r'p(r')d"^r' is the electric dipole moment. The LIV terms in (k^-)^'' break 
the radial symmetry giving a non-Coulomb behavior to the static solutions. Despite the large number of terms 
in these solutions, we verify that the electric field preserves the and the decaying behaviors for the 
monopole and dipole moments, respectively, as occurs in the pure Maxwell electrodynamics. Obviously, it is a 
consequence of the dimensionless character of the LIV coefficients. 

B. The Green's function for the magnetic field 

Now, we search for an explicit solution for the magnetic field. The Green's function for the magnetic field 
equation of motion (jlip is written as 



(1 - n)5ai - [ne-f ) + {H.^y dad, Gib (r - r') = 5ab5\r - r') 



(25) 



Using the Fourier representation and having much care in the tensor inversion procedure, we obtain in the 
momentum space 



Gab (p) = — \ 



(1 + n) Sab + iKc-)ab - ('«e-)c6 



PaPc 



(26) 



Performing the inverse Fourier transformation, we attain the following expression: 



Gabir - r') = I (1 + n) Sab + %^ + ^^"^^ ' ^'^^^^ ' 



4tt |r — r\ 

with which the magnetic field is then written as 



2|r- r'l 



(r) = - / d'r' G,, (r ~ r') (V x J (r'))'' 



(27) 



(28) 
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It leads to the explicit solution: 



(r) 



.3^, (Vxj(rOr 



(V X J (r'))' 



(r-r'),(r-r'r 



After a certain algebraic effort, a dipolar expansion for the magnetic field is achieved as well, yielding 



(r) 



1 

47r 



(1 + n) 



TO 



3 (r • m) 



3 m* 



15 (r • m) 



(29) 



(30) 



where we have considered a localized and divergenceless current density distribution J, and m =i J r' x J(r')d^r' 
is the magnetic dipole moment. In Eq. ([50]) the first term inside the parentheses is the usual Maxwell contri- 
bution, just corrected by the (1 + n) factor. The terms that are proportional to the LIV coefRcients, {ke-Y^ ■, 
ascribe to the solution a directional dependence or anisotropic character. In principle, such a directional depen- 
dence could be used to impose an upper bound on the LIV parameters. In Sec. Ill, we show that the attained 
bound is not as restrictive as desired. 



III. FINAL REMARKS 



We should now compare these parity-even stationary solutions with the parity-odd ones derived in Ref. [11] . 
At the stationary regime, the main difference is that now the electric and magnetic sectors are not coupled by 
the LIV tensor anymore. In the parity-odd case, a stationary current is able to produce an electric field as much 
as a static charge can generate a magnetic field. As such an interconnection does not appear in the present case, 
the manifestation of pure LIV electromagnetic effects (aside from Maxwell ones) , as the production of magnetic 
field by a static point-like charge (see Ref. ll|), are absent. Now, the LIV effects appear as small corrections for 
the usual Maxwell's electric and magnetic fields. Yet, the LIV effects can still be identified by means of suitable 
devices, as it is discussed below. Apart from this difference, the solutions of the parity-even and parity-odd 
sectors possess some similarities. Indeed, the electric field for a point-like charge (in both sectors) exhibits an 
asymptotic behavior as , while a stationary current provides a magnetic field whose dipolar expansion is 
proportional to mr~^ . This is ascribed to the dimensionless character of the tensor W . 

The attained magnetic field solution docs not lead to good upper bounds on the magnitude of the parameters 
n and {lie-)n) when we take as reference the Earth's magnetic field. In fact, proceeding in a similar way as in Ref. 

we assert that the LIV tensor must not imply a magnetic field contribution larger than 10~^G (otherwise 
it would be detected). From Eq. ([HO)) , we observe that the LIV terms are always proportional to (m/r^). 
Assuming that to represents the Earth's magnetic dipole, and the Earth's radius, it holds the following 
ratio m/R^ = 0.3 G (see Ref. Q). This procedure, however, implies a non-restrictive bound: n < 10 . 

A much better bound for the parameters (we-)*"' can be attained from the expression for the scalar potential. 
The idea is to evaluate the scalar potential generated by a charged sphere in different outer points located at the 
same distance from the center of the sphere, observing the difference of potential induced by the non-Coulomb 
LIV term. The starting point is the expression for the potential generated by a conducting sphere of radius R 
and charge q (uniformly distributed over its surface), which can be achieved by replacing the charge density for 
a sphere, p(r') = q5{r' — R) / {in R'^) , in Eq. Using Fourier integrations (see the appendix), the potential 

is (for r > R) 



Ao (r) 



47r r 2 



ran r 



(31) 



We see that the term in (Ke-)'"* breaks the radial symmetry of the potential, implying potential variations along 
a circular path around the center. We now expand the term (kc-)"'' ''a^'b at the form 

(ke-T'' ran ^ (Ke^)^^ {rif - (rsf + (ke^f^ {r2f ~ {r^) 



+ 2 (ke-)^^ rir2 + 2 (ke-)" rira + 2 {ne-Y'" r2r'3. 



n23 
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where we have used the traceless matrix 



(^e-)ii ('*e-)i2 
(^e-) = I (Ke-)i2 ('*e-)22 
('*e-)l3 ('*e-)23 



('«e-)i3 
('^e— )ii ~ ('^e-)22 



(32) 



Then, we can conceive of an experiment to measure the electrostatic potential generated by a 1 C charged 
sphere of radius R (maintained in vacuum) in two distinct outer points, A and B, located at the a circle of 
radius r > R on the x ~ y plane. We consider the points A and B symmetrically disposed in relation to the 
y-axis at the positions: A — r(cos (j), sin (p, 0), B — r{— cos </>, sin </>, 0). Then, the difference of potential between 
these points is simply 



AAo = Ao{A) ~ Ao{B) = -f {K,.y^ sin 20 (r^_R^ 
for r > R. For (p — tt/A and a 1 C charge, such difference of potentials equal to 



A^o = 9 X 10^ (Ke-) 



12 



(33) 



(34) 



For attaining the best bound, we should consider the maximum value of Eq. (|34p . So, it must be evaluated 
at the point r = Ry/3, in which the expression (r^ — i?^) r^"^ has a maximum. For a charged sphere of 
unitary radius {R = Im) , we obtain AAq = 3.46 x 10^ (ke-)"'^^ V. Given the existence of sensitive methods for 
measurement of the potential able to detect slight variations of 1 part in 10^° V, we can infer that the voltage 
difference of Eq. ((34)) cannot be larger than lO^^*' V, that is, 3.46 x 10^ (ke-)^^ < lO^^*'. This condition leads to 

z and x~ z, this upper limit holds equivalently 



{Ke-) < 2.9 X 10 Choosing pairs of points on the planes y ^ 



for (ks-)'^^ and (^e-)^^. 



This device can also be used to set up an upper bound on the diagonal components (kc- 
(ke-)^^) we take the points A and B in the positions: A = r(l, 0, 0), B 
between these points is 



For constraining 
r(0, 0, 1). The difference of potential 



AAo ^ Ao{A) - Ao{B) = ± (Ke-)'' 

An 



(35) 



which leads to the same bound obtained for the non-diagonal components: (ke-)^^ < 2.9 x lO"'^'^. Choosing two 
points on the y — z plane, A — r(0, 1,0), B — r(0, 0, 1), this bound can be stated to (kc-)^^ . Thus, we conclude 
that by means of this experiment it is possible to establish an upper bound as stringent as 



(Ke-Y^ < 2.9 X 10 



-20 



for the five non-isotropic components of the traceless matrix (Kg-)'"' . This is a nice bound for an Earth-based 
experiment, as good as the best bounds stated from astrophysical data analysis of UHECRs [16]. 

As the isotropic component n does not break the spherical symmetry of the potential, this kind of experiment 
does not provide any way for bounding it. This component induces a slight screening on the Coulomb potential 
that may be interpreted as a charge screening. An experiment able to constrain n could be based on a charge 
or potential screening measurement. In this case, the major difficult is that the tiny LIV effect is disguised by 
the dominant Maxwell's behavior, avoiding its isolation. So, the LIV effect stays limited by the experimental 
imprecisions of the device. A two-sided bound was recently stated for this coefficient in the context of quantum 
electrodynamics decay processes modified by this LIV parameter 



Finally, we should note that this work completes the calculation of the stationary solutions of Maxwell's 
electromagnetism modified by the non-birefringent elements of the abelian CPT-even and LIV sector of the 
standard model extension, a task initiated in Refs. lol. 11|. 
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APPENDIX A: EVALUATION OF THE SCALAR POTENTIAL GENERATED BY A CHARGED 

SPHERE 

By starting from Eq. (jl9p . the scalar potential is rewritten as 



where 



..w.±{a-..)/«v^-2fc)!!.4 



r — r 



(Al) 



(A2) 



and using (Ke_)j,j = tr (ke-) — 0, we show that 



,~ sab (r - r )a{r ~ r )b _ _ .ab ^ 

l^e-j ; -3 — l'*e-j 

|r — r'l L I- - I J 

Knowing that the charge density for a charged sphere of radius R is p{r') = q5 (r' — R) /{AttR'^), its Fourier 
transform is 



sin (pR) 
pR ' 



9- 



with p = |p|. For evaluating the integral lb, we use the Fourier representation of the Coulomb potential 



r — r 



47r 



— ip- (r — r/) 



(27r)^ p2" 



(A3) 



(A4) 



By substituting in the Eq. (|A2p and using (jASp we obtain, after some algebraic manipulations, the following 
expression for Ii,: 



Rj {2TTf VP 



^^'V (27r)^ VP^ 



p^ J 



d^n ( 1 



{2Tif VP^ 



COS 



Solving these integrals, we obtain 



,(r-r')f, , g qR"^ 



—p{v) = -rb 
r — r' r 



3r3 



Thus, finally, we obtain the scalar potential generated by the charged sphere of radius R for (r > R) 



ra^ft (r^ - R? 



(A5) 



(A6) 



(A7) 
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